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Abstract 

A direct calculation of the mixed-action parameter A mix with valence overlap fermions on a domain- 
wall fermion sea is presented. The calculation is performed on four ensembles of the 2+1-flavor 
domain-wall gauge configurations: 24 3 x 64 (ami = 0.005, a = 0.114 fm) and 32 3 x 64 (ami = 
0.004,0.006,0.008, a = 0.085 fm). For pion masses close to 300 MeV we find A mix = 0.030(6) GeV 4 
at a = 0.114fm and A mlx = 0.033(12) GeV 4 at a = 0.085 fm. The results are quite independent 
of the lattice spacing and they are significantly smaller than the results for valence domain-wall 
fermions on Asqtad sea or those of valence overlap fermions on clover sea. Combining the results 
extracted from these two ensembles, we get A m « = 0.030(6) (5) GeV 4 , where the first error is 
statistical and the second is the systematic error associated with the fitting method. 

PACS numbers: 12.39.Fe, 11.30Rd, 12.38.Gc 



I. INTRODUCTION 

Mixed action approaches have been studied by several 
groups, such as Domain-wall fermion (DWF) valence on 
Asqtad fermion sea [1], overlap valence on DWF sea [2], 
overlap valence on clover sea [3], and overlap valence on 
twisted-mass fermion sea [4] . In view of the fact that it is 
numerically intensive to simulate chiral fermions (DWF 
or overlap), it is deemed practical to use the cheaper 
fermion formulation for generating gauge configurations 
and the more expensive fermion discretization for the va- 
lence as an expedient approach toward dynamical QCD 
simulations with chiral fermions. Many current alge- 
bra relations depend only on the chiral properties of the 
valence sector. The mixed action theory with different 
fermions for the valence and the sea is a generalization of 
the partially quenched theory with different sea and va- 
lence quark masses. The mixed action partially quenched 
chiral perturbation theory (MAPQxPT) has been devel- 
oped for Ginsparg- Wilson fermions on Wilson sea [5] and 
staggered sea [6], and has been worked out for many 
hadronic quantities to next-to-leading order (NLO), such 
as pseudoscalar masses and decay constants [5-7] , isovec- 
tor scalar ao correlator [8-11], heavy-light decay con- 
stants [12], and baryon masses [11, 13]. 

In the mixed action chiral perturbation theory with 
chiral valence fermions, it is shown [5] that to NLO there 
is no 0(a 2 ) correction to the valence- valence meson mass 
due to the chiral symmetry of the valence fermion. Fur- 
thermore, both the chiral Lagrangian and the chiral ex- 
trapolation formulas for hadron properties to the one- 
loop level (except ^-dependent quantities) are indepen- 
dent of the sea fermion formulation [14]. The LO mixed- 
action chiral Lagrangian involves only one more term 



with 0(a 2 ) discretization dependence which is charac- 
terized by a low energy constant A m i X . The LO pseu- 
doscalar meson masses for overlap valence and DWF sea 
are given as 

m 2 vv , = B ov (m v +m v >), 

m 2 vs = B ov m v + B dw (m s + m res ) + a 2 A mixi 
m 2 ss , = B dw {m s + 77V + 2m res ), (1) 

where m vv > / m ss > is the mass of the pseudoscalar meson 
made up of valence/sea quark and antiquark. m vs is the 
mass of the mixed valence and sea pseudoscalar meson. 
Up to numerical accuracy, there is no residual mass for 
the valence overlap fermion. The DWF sea has a residual 
mass m res which vanishes as Ls — > oo. A m j a enters 
in the mixed meson mass m vs as an 0(a 2 ) error which 
vanishes in the continuum limit. We should note that, 
unlike the partially quenched case, even when the quark 
masses in the valence and sea match, the unitarity is still 
violated due to the use of mixed actions. The degree 
of unitarity violation at finite lattice spacing depends on 
the size of A mix . 

A m i X has been calculated for DWF valence and Asq- 
tad fermion sea which gives A m « = 0.249(6) GeV 4 [15] 
at a = 0.125 fm, 0.211(16) GeV 4 at a = 0.12 fm [10], 
and 0.173(36) GeV 4 at a = 0.09fm [10]. It is also cal- 
culated for overlap valence and clover sea which yields 
A mix = 0.35(14) /0.55(23) GeV 4 at m w = 190/300 MeV 
and a = 0.09 fm [16]. This means that for a valence pion 
of 300 MeV, A mix produces, at a = 0.12 fm, a shift of 
- 110 - 240 MeV and, at a = 0.09 fm, a shift of ~ 
55 — 153 MeV for these cases, which are substantial por- 
tions of the valence pion mass. 

We have used valence overlap fermions on the 2+1- 
flavor DWF sea to study hadron spectroscopy [17-20]. 
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In this work, we calculate A m ix for such a mixed ac- 
tion approach which is needed for chiral extrapolation in 
MAPQxPT. 



for all the masses in each ensemble. For the 24 3 x 64 
ensemble the details are presented in Table II. For one of 
the sea quark mass in the 32 3 x 64 ensemble, the results 
are presented in Table III. 



II. CALCULATION DETAILS 

The A m i x parameter is calculated on four ensembles of 
the 2 + 1-flavor domain wall fermion gauge configurations 
[21, 22]. Two different lattice spacings were used to study 
the dependence on the cutoff. In addition, we also used 
multiple sea masses, for the 32 3 x 64 lattices, to study 
the sea quark mass dependence of A mix . Details of the 
ensembles are listed in Table I. 



Lattice Size 


a -1 (GeV) 


ami 


m ss (MeV) 


24 3 x 64 


1.73(3) 


0.005 


329(1) 


32 3 x 64 


2.32(3) 


0.004 


298(1) 


32 3 x 64 


2.32(3) 


0.006 


350(2) 


32 3 x 64 


2.32(3) 


0.008 


399(1) 



TABLE I: Details of the DWF ensembles used in this work. 

There arc different parameterization schemes [10, 15] 
one can use to relate A m i X to the quark and pseudoscalar 
masses from Eq. (1). In this paper we choose to param- 
eterize A m i x as 

m ls - \ m2 ss = B o V m v + a 2 A mix . (2) 

This is similar to the parameterization used in [10]. The 
quantity Sm 2 = m 2 s — \rn 2 ss , has a linear behavior in m v 
and can be calculated directly by computing only the 
pseudoscalar masses m vs , and m ss . We see that in the 
regime where Eq. (2) is valid, A« is equivalent to dm 2 
in the limit m v — > 0. 

In this work we used an overlap operator with a HYP- 
smeared kernel. This was shown to have better numer- 
ical properties [17]. We calculated the masses m vs and 
m ss using 50 configurations for each ensemble. Recall 
from Section I that m ss requires the propagators for the 
sea quark mass which were computed with DWF. m vs 
needs the propagators for both overlap and DWF. The 
DWF propagators were made available by LHPC. The 
overlap propagators were computed using a polynomial 
approximation [23] to the matrix sign function. They 
were used to compute 14 — 16 values of m vs . A multi- 
shifted version of the conjugate gradient algorithm [24] 
was implemented for the overlap propagator calculation 
to compute all masses at once. To accelerate the inver- 
sions for the overlap propagators we employed a deflation 
technique which has been seen to speed up the calcula- 
tion significantly [17]. 

We fitted correlators using single state exponential fits 
to extract the pion masses. The fitting windows were 
adjusted to get a reasonable y^/dof and are the same 



TABLE II: Pion mass fitting details for the 24 3 x 64 lattice. 
m ss ~ 329 MeV. 



am v 


m„ s (MeV) 


xls/dof 


m„„(MeV) 


xlv/dof 


0.0014 


274(3) 


1.0 


122(4) 


1.8 


0.0027 


278(2) 


1.1 


154(2) 


1.4 


0.0046 


287(2) 


1.2 


188(2) 


1.6 


0.0081 


305(2) 


1.2 


242(2) 


1.6 


0.0102 


315(2) 


1.2 


270(2) 


1.4 


0.0135 


331(2) 


1.1 


308(2) 


1.3 


0.0153 


339(1) 


1.1 


327(2) 


1.3 


0.0160 


342(1) 


1.1 


334(2) 


1.3 


0.0172 


347(1) 


1.1 


346(2) 


1.3 


0.0243 


378(1) 


1.2 


409(1) 


1.6 


0.0290 


397(1) 


1.2 


445(1) 


1.7 


0.0365 


426(1) 


1.2 


498(1) 


1.5 


0.0434 


451(1) 


1.2 


542(1) 


1.3 


0.0489 


471(1) 


1.2 


576(1) 


1.2 


0.0670 


531(1) 


1.0 


677(1) 


1.3 


0.0710 


543(1) 


1.0 


698(1) 


1.4 



TABLE III: Pion mass fitting details for the 32 3 x 64 lattice 
for ami = 0.004. m ss ~ 298 MeV. 



am v 


m„ s (MeV) 


xls/dof 


m„„(MeV) 


xlv/dof 


0.0007 


237(3) 


1.0 


126(10) 


1.5 


0.0015 


238(2) 


1.0 


136(7) 


1.2 


0.0025 


246(2) 


1.0 


160(4) 


1.2 


0.0035 


254(2) 


1.0 


184(3) 


1.2 


0.0046 


264(2) 


1.0 


209(3) 


1.2 


0.0059 


274(2) 


1.0 


235(3) 


1.2 


0.0068 


282(2) 


1.0 


253(3) 


1.3 


0.0076 


289(2) 


1.0 


268(3) 


1.4 


0.0089 


299(2) 


1.0 


288(3) 


1.6 


0.0112 


317(2) 


1.1 


324(4) 


1.0 


0.0129 


329(2) 


1.1 


337(5) 


1.0 


0.0152 


344(2) 


1.2 


365(5) 


1.3 


0.0180 


362(2) 


1.3 


396(6) 


1.5 


0.0240 


399(3) 


1.6 


454(9) 


1.9 



III. FITTING STRATEGIES AND RESULTS 

Extracting A mix requires fits to squares of multiple 
pion propagators with different valence quark masses; the 
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FIG. 1: (am m ) 2 /2( am v ) as a function of am v for the 24 3 x 64 ensemble (left) and the three 32 3 x 64 ensembles (right). 



resulting masses will be correlated since they all come 
from a single ensemble. To compute the covariance ma- 
trix we would need to perform an augmented x 2 ~fit m ~ 
volving the pion propagators for all valence quark masses 
simultaneously. This is not numerically stable due to the 
large number of parameters in the model. We thus have 
to use an alternative procedure to take into account these 
correlations. To gauge the systematic errors introduced 
by our choice of fitting method we will use multiple fitting 
procedures. In this section we will describe three differ- 
ent fitting strategics. These methods differ in the way 
we account for correlations among the different valence 
masses. 

In method I we follow the standard jackknife philoso- 
phy by defining a A m i X estimator directly in terms of the 
raw pion propagators; the error is then determined by the 
variance over the jackknife ensemble. In method II we use 
the jackknife procedure to estimate the Sm 2 covariance 
matrix and then compute A m i X using a standard corre- 
lated fit [25]. These two procedures are used for both the 
24 3 x 64 and 32 3 x 64 ensembles. However, since there 
are three sea quark masses for the 32 3 x 64 volume, we 
can perform an uncorrelated fit using Sm 2 computed on 
the three independent ensembles — this is our method III. 
From the different methods we can obtain a systematic 
uncertainty for our results. 

In each method we performed a binned jackknife anal- 



ysis of Sm 



1 2 

2 III e o 



Because we performed four 



inversions with four different point sources per configu- 
ration we binned in units of four and eight. This allowed 
us to drop one or two whole configurations in each jack- 
knife subsample. We find the uncertainties in both cases 
to be of comparable size, indicating that autocorrelation 
is negligible. 

The first step in our fitting is to determine a range 
of quark masses where the tree-level relation between 
m 2 v and m v holds so that we can use Eq. (1). Be- 
low this range, one expects to see chiral logs including 
partially quenched logs and other non-linear m v depen- 
dence from the NLO in xPT. Above this range, tree-level 



%PT is not expected to be valid. We do this by plotting 
(am vv ) 2 /2{am v ) as a function of am v as shown in Fig. 1. 
We choose the fitting range in the region where the ratio 
(am vv ) 2 /2(am v ) is fairly flat; these ranges are tabulated 
in Table IV and are used in all three fitting methods. We 
note that in the range we are fitting, m^L > 4 for both 
m vs and m vv so that the volume dependence is expected 
to be small. 



A. Method I: weighted averaging 



In this method we used a weighted linear fit to extract 
the value of A mix for each bin. The weights, <J 2 m2l are 
given by 



(3) 



where <J m ^ s , and u mss are the uncertainties of the mixed 
and DWF pion masses respectively. Eq. (3) was derived 
using the standard error propagation formula neglecting 
correlation between m vs and m ss . The cross-correlations 
will be accounted for by the external jackknife procedure. 
By using weighted fitting, less importance is given to data 
points with larger uncertainties. After fitting each bin we 
then have a jackknife ensemble {A m i X }. We use square 
brackets, [ ] , to indicate a particular jackknife sample and 
angle brackets, ( ), to denote jackknife averages. For the 



Lattice 


ami 


am v fit range 


24 3 64 


0.005 


0.0243 - 0.0489 


32 3 64 


0.004 


0.0112 - 0.0240 


32 3 64 


0.006 


0.0112 - 0.0240 


32 3 64 


0.008 


0.0112 - 0.0240 



TABLE IV: Range of quark masses, am v , used in the fitting 
procedures to extract A mix via Eq. (1). 
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FIG. 2: Extracting A m i x from a linear extrapolation of 8m 2 for the 24 3 x 64 ensemble (left) and the three 32 3 x 64 ensembles 
(right). For the 32 3 x 64 plot the inset shows the intercept of the fit and the error bars of the extracted values of A m i X . 



case of A. m i Xl its jackknife average and the corresponding 
uncertainty is given by 



k=l 



*A mi „ = ^(N-l)((A 2 mix )-(A mix )^). 



(4) 



We minimize the correlated x 2_ f un ction, 

x 2 = J2 «M<) - MO C *J 

x ( ({5m 2 ),) (Sm%) , 



with the covariance matrix, Cj given by the jackknife 
estimate 



The extracted values of A m ,- X , using this method, are 
presented in Table V. We also list the corresponding val- 
ues of B ov . Figure 2 shows a 2 8m 2 as a function of am v 
and the corresponding linear fit. 



B. Method II: covariance matrix 



In our second method we perform a correlated valence- 
quark mass fit to the function 5m 2 = B ov m v + a 2 A mix . 
The central values are taken to be the jackknife averages, 



(N-l) 
N 



N 



£([M<]*-<M<» 



fe=i 

x (W] k -(Mi» 



(6) 



The subscripts i and j index the valence-quark mass. The 
factor (N — l)/N differs from the usual definition of the 
covariance matrix to account for the correlation of the 
jackknife samples [25]. The fit uncertainties are obtained 
by constructing the standard Hessian matrix. Results of 
this method are tabulated in Table V. 



(Mi>=»;£[Mi] fc - (5) 

k=l 



Lattice 


ami 


A m „(GeV 4 ) 


B ov {GeV) 




I 


II 


I 


II 


24 3 64 


0.005 


0.032(6) 


0.028(5) 


1.85(2) 


1.88(2) 


32 3 64 


0.004 


0.040(14) 


0.025(9) 


1.88(10) 


1.95(6) 


32 3 64 


0.006 


0.054(9) 


0.050(8) 


1.81(5) 


1.81(5) 


32 3 64 


0.008 


0.059(13) 


0.063(13) 


1.74(6) 


1.73(6) 



TABLE V: Extracted values of A m i X and B ov using fitting 
methods I and II. 



C. Method III: uncorrelated fitting 

The methods described in the two previous sections 
closely parallel the method performed by [10] in which 
for each lattice ensemble a string of partially quenched 
meson masses are calculated; from them we extract A mix . 
In this section we perform a fit on the three 32 3 x 64 
ensembles based on the value of Sm 2 measured at the 
point where the valence and sea quark masses match [26] . 
There are no cross-correlations among the masses since 
the ensembles are independent. 

The value of 5m 2 is computed for the pion mass that 
most closely satisfies the condition m ss s» m vv . Because 
it is not easy to match a priori the pseudoscalar masses 
and it is expensive to regenerate overlap propagators with 
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FIG. 3: Determining A m i! by performing the linear extrap- 
olation as in Fig. 2. Each point corresponds to one of the 



32 x 64 lattices where m a 



A magnified view of the 



extrapolation in the neighborhood of the chiral limit, along 
with the extrapolated value, is shown in the inset. 
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FIG. 4: A m i X for a — 0.085 fm. The empty symbols indicate 
the results of method I and the full symbols are from method 
II. The continuous line and the shaded region are the results 
of the chiral extrapolation in method III. 



different masses, we performed interpolation among the 
existing data points to obtain a better approximation of 
where m ss and m vv match. 

We do an uncorrelated y^-ftt with the error bars ob- 
tained by a jackknife procedure. Figure 3 shows the 
data points and the fit results. For this case we find 
A mlx = 0.042(24) GeV 4 , and B ov = 1.82(16) GeV. 



IV. DISCUSSION 

As a first step, we look at the lattice spacing depen- 
dence. We compare the two ensembles with m ss close 
to 300 MeV. These correspond to am; = 0.005 and 
ami = 0.004 with a = 0.114 fm and a = 0.085 fm re- 
spectively. We average the central values and the errors 
from the two fitting methods, separately for each lattice 
spacing. We find at a = 0.114 fm A rnlx = 0.030(6) GeV 4 
and for a = 0.085 fm A mlx = 0.033(12) GeV 4 . We see 
that the lattice spacing dependence is smaller than our 
errors; this indicates that A m i X is capturing the dominant 
lattice artifact for the parameters used in this study. 

To discuss the sea quark dependence of A mix we plot 
in Fig. 4 the results of methods I, II, and III for the 
ensembles with a = 0.085 fm. We note that the results 
are consistent within two sigma. This is consistent with 
LO MAPQxPT in Eq. (1) where A mix is a low energy 
constant, independent of the valence and sea masses. 

We now combine the A m j X values extracted from the 
two lightest sea mass ensembles to produce our final re- 
sult. We use these ensembles because the LO MAPQxPT 
is expected to describe the data better at lower sea quark 
masses. For each of the fitting methods we combine the 
results of the a = 0.114 fm and a — 0.085 fm ensembles. 
Since the two ensembles are statistically independent, it 
is straightforward to combine these results: for method I 



we get A m i X = 0.033(6) GeV 4 and for method II we get 
= 0.027(5) GeV 4 . We can now use the values de- 
termined using these two methods to estimate the sys- 
tematic fitting errors. We quote the final result with two 
uncertainties, the first statistical and the second associ- 
ated with fitting systematics. The central value and the 
statistical error are taken to be the average of the results 
from methods I and II. The systematic error is the stan- 
dard deviation of the results from these two methods. 
We get A mix = 0.030(6)(5) GeV 4 . 

Table VI lists calculated values of A m i X for pion masses 
close to 300 MeV using different mixed actions. We no- 
tice that our values of A m j a are significantly smaller 
than overlap on clover or DWF on Asqtad for compa- 
rable lattice spacings and pion masses. For both m vv 
and m ss at 300 MeV, our calculated A m i X will shift the 
pion mass up by 10 and 16 MeV for the 32 3 x 64 lattice 
at a = 0.085 fm and 24 3 x 64 lattice at a = 0.114 fm, 
respectively. These are substantially smaller than the 
corresponding 55 - 153 MeV and 110 - 240 MeV shifts 
that we mentioned in Sec. I. The value of A m i X is sig- 
nificantly smaller in our case probably because the sea 
and valence fermion actions are similar. Both of them 
are approximations of the matrix sign function, but with 
different kernels. 

The previous calculation [17] of A m i X , for overlap on 
the DWF sea, has roughly the same value but the sign is 
negative. That calculation measured A mix by examining 
the states that wrap around the time boundary. This 
indirect method is less reliable and the errors are large. 



V. CONCLUSION 

We calculated the additive mixed action pseudoscalar 
meson mass parameter, A mix , for the case of valence 
overlap fermions on a DWF sea. A ralI is significant be- 
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Mixed Action 


Ref. 


a fm 


A mix (GeV 4 ) 


DWF on staggered 


[15] 


0.125 


0.249(6) 


DWF on staggered 


[10] 


0.12 


0.211(16) 


DWF on staggered 


[10] 


0.09 


0.173(36) 


overlap on clover 


[16] 


0.09 


0.55(23) 


overlap on DWF 


this work 


0.114 


0.030(6) 


overlap on DWF 


this work 


0.085 


0.033(12) 



TABLE VI: A m i X values for DWF valence quarks on stag- 
gered sea quarks for pion mass at 300 MeV. 



dence of A m ix at a = 0.085 fm and we find that they 
agree within two sigma. Combining the results extracted 
from the ensembles with the lightest sea quarks, we get 
A m i X = 0.030(6)(5) GeV 4 , where the first error is sta- 
tistical and the second is the systematic error associated 
with the fitting method. 

When compared to previous mixed action studies, 
DWF on staggered or overlap on clover, the values of 
Amix of overlap on DWF are almost an order of magni- 
tude smaller. This is most likely due to the fact that the 
sea and valence fermions are similar. 



cause it enters into mixed action partially quenched chiral 
perturbation theory (MAPQxPT) for chiral extrapola- 
tion of many low energy observables. 

Two different lattice spacings were used to examine 
the cut-off behavior. For a pion mass close to 300 MeV 
we find A mlx = 0.030(6) GeV 4 at a = 0.114fm and 
A mlx = 0.033(12) GeV 4 at a = 0.085 fm. They are the 
same within errors. Our calculated A m i X will shift the 
pion mass up by 10 and 16 MeV for the 32 3 x 64 lattice 
at a = 0.085 fm and 24 3 x 64 lattice at a = 0.114 fm, 
respectively. We studied the sea quark mass depen- 
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